Crater Property in Two-Particle Bound States: When and Why 
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Crater has shown that, for two particles (with masses mi and 7712) in a Coulombic bound state, 
the charge distribution is equal to the sum of the two charge distributions obtained by taking 
mi ^ 00 and m2 —> cx3 respectively, while keeping the same Coulombic potential. We provide a 
simple scaling criterion to determine whether an arbitrary Hamiltonian possesses this property. In 
particular we show that, for a Coulombic system, fine structure corrections preserve this Crater 
property while two-particle relativistic corrections and/or hyperfine corrections may destroy it. 



Recently, in an interesting paper in this journal M, Crater discussed an unusual feature of charge densities for 
two-particle Coulombic bound states. Let p{R] mi, 7712) be the charge density of a two-particle bound state in a given 
potential V{r) = V{ri —72) in the center-of-mass coordinate system. Then Crater observed that, for a Coulombic 
potential the charge density satisfies the following relation: 



p{R;mi,m2) 



lim p{R; mi, 1712) + lim p{R;m,i,m,2), 

7n2 — ^cxD nil — 'OO 



(1) 



or, in Crater's own words, one can picture p{R', mi, 7712) "as equivalent to that produced by a particle of mass ?ni and 
charge ei, bound to a fixed center with charge 62 plus that produced by a particle of mass m,2 and charge 62, bound 
to a fixed center with charge ei," with the fixed center being the center of mass. Here and after this property will be 
referred to as the Crater property. Crater has shown in Ref. [Q that the Crater property holds for Coulombic potentials 
but not for generic potentials. In the real world, however, Coulombic potentials are often corrected by perturbations 
like fine /hyperfine structures and relativistic effects. It would be of interest to know what kind of corrections can be 
added to the Coulombic potential without destroying the Crater property. More generally, we would like to have a 
criterion to determine whether a given potential has the Crater property without explicitly solving the Schrodingcr 
equation. The purpose of this paper is to provide answers to these questions. 

Let us consider an eigenfunction ^[r;m,i,m2) to the Schrodingcr equation (in units h = 1): 



H{mi, m2)i^{r; mi, 7712) 






V{r;mi,m2) 



ijj{r;mi,m2) = E{mi, m2)ip{r; mi, m2), 



with /_* = mim2/M and M = mi + m2- The charge density operator p{R) is defined as |I| 

p{R) = eiS'^iR - n) + e2S^{R - 72), 
and the charge density p{R; mi, 77*2) is its expectation value, which can easily be shown to be 
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Since the eigenfunction 7/;, satisfying the normalization condition Jd'^rji/'p = 1, carries scaling dimension 
[length] ~'^/^ = [momentum] ■^/^ = [mass]'^/^ (with units Tt. = c = 1), it is always possible to rewrite ■(/'(r; 7771, 7772) 
as 



7/)(r; 7771, 7772) = ^^/^7/J(r; 7771, 7772) 



(5) 



where the rescaled eigenfunction ip{r;mi,m,2) is a dimensionless function. Now consider the case when 7/; (r; 7711,7712) 
has the following form: 



7/^(7; 7771, 7772) = ipilir), 



(6) 



which states that the dependences of ip on the location r and the masses toi_2 always come through the combination 
Hr = TOi?n2r/(mi +1712). Then ?/'(r; toi, 7712) — /i'^'^?A(/i?') and the charge density p{R;mi,m2) in Eq. (]^ becomes 
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lim p(i?; TOi, 7712) 

m2 — ^00 



lim p(R;mi,m2) 

mi — >oo 



(7) 



which is exactly the expression for the Crater property. In other words, the charge density exhibits the Crater property 
whenever the eigenfunction can be written as p,^''^ip{p,r). 

It is easy to translate the above scaling condition on the eigenfunction to a corresponding scaling condition on 
the Hamiltonian. Since the Hamiltonian H(mi,m2) carries scaling dimension [mass]^, it can always be rewritten as 
/iiJ (toi, 7712), where 7f is a dimensionless function of ?7ii.2, as well as the relative coordinates r and the canonical 
momenta —id /dr. It is straightforward to see that V' is a function of solely p,r if and only if 



H{m.i,m2) —Hifir. 



d 



difir) 



+ Vo, 



(8) 



such that, up to an additive constant, the masses enter the Hamiltonian only through combinations /ir and 
—i{d/d{fj,r)). The dimensionless constant Vb, which may have arbitrary dependences on ?ni_2, niay shift the eigenvalues 
but does not affect the eigenfunctions. 

We have shown that the charge density of an eigenfunction exhibits the Crater property if and only if the Hamiltonian 
can be written as 



H ^ p 



n[ pr,- 
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which will be referred to as the scaling criterion. With this criterion one can easily determine if a particular potential 
exhibits the Crater property. For spinless Schrodinger systems, the kinetic term always satisfies the scaling criterion. 
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On the other hand, for analytic V^(r)'s one can expand them in Laurent series and the scaling criterion is satisfied if 
and only if 
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+ Vo., 



(11) 



where ak are mass independent coefficients. Of special interest is the case where the only non-vanishing Ofc are those 
with fc = — 1 and —2: 
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which describes a Coulombic potential in three dimensions with a_i — 6162 and a-2 ~ £{£+ 1)/2, i.e., the case studied 
in Ref. &. Another interesting case is the "two-dimensional Coulombic potential", i.e., the logarithmic potential: 



V{r) = 6162 ln(r/ro) = £162 [In(^r) - ln(/iro)] 



(13) 



which can be seen to satisfy the scaling criterion by identifying the second term as Vq. 

As pointed out in Ref. |l||, the Crater property is not a feature of potentials of generic r and mass dependences. 
Crater illustrated this point by studying the eigenfunctions of a simple harmonic potential and showed explicitly that 
p{R;mi,m2), given by Eq. (0), is not the sum of limmj^oo p(i?; 7711,7712) and lim„jj_^oo p(^; "1ij"^2)j J-e., the simple 
harmonic potential does not exhibit the Crater property. On the other hand, we can reproduce the same conclusion 
by just studying the scaling behavior of the simple harmonic potential: 



V{r) = -M^V, 



(14) 



which cannot be recast in a form conforming to criterion (||). As a result, the Crater property is not exhibited in 
simple harmonic potentials. 

It is of interest to note that, for any potential V{r) — ^V{iir) satisfying the scaling criterion, including the fine 
structure corrections does not destroy the Crater property. 



1 ldV{r) 1 1 dV{^Jir) ^ 

^^c^ r ar & /ir d[fir) 

1 d^Vjr) ^ 1 d^Vjr) 

where Hmv, Hso, a-nd Hjj stand for the relativistic mass variation term, the spin-orbit coupling term and the Darwin 
term, respectively. This may look miraculous but is actually nothing but a consequence of the fact that all these fine 
structure corrections come from the leading order expansion in the fine structure constant of the one-particle Dirac 
Hamiltonian with the same potential V{r). 



d 

H = -ia ■ — +l3^i + V{r) = n 
or 



.a.^-H/3-Hy(H 



(18) 



where a and /3 are the Dirac matrices. Since this one-particle Dirac Hamiltonian also satisfies the scaling criterion, 
the Crater property is preserved. 

However, it is important to bear in mind that it is an approximation to describe a two-particle bound state by a 
one-particle Schrodinger or Dirac equation. Take, for example, the hyperfine correction, which for a Coulombic bound 
state is 



51526162 3 



3toiTO2 



M 



M 



51526162 



H,f = ^i^^^^^S, ■ 52 S^r) = ^ ^ £i£frin5, . S, 6^^,r), (19) 



1 92 1 92 


1 92 ^ Q2 


H 9 


— 1 


2toi drf 2to2 ^rj 


2M(9i?2 ' 2/i9r2' 



and the scaling criterion is violated by the outstanding factor of /i/Af , where M ~ mi + rn2 is the total mass. 
Violations of the scaling criterion may also be due to two-particle relativistic effects. In the non-relativistic theory 
a two-particle problem can always be reduced to an effective one-particle problem in the relative coordinates, which 
decouple with the center-of-mass coordinates in the kinetic term: 

I 1 «2 

(20) 

where R is the center of mass position miri + m2r2, and r is the relative position ri — r2. For a relativistic theory 
in general no such decomposition is possible, and the description of a two-particle problem by a one-particle equation 
is a good approximation only when one particle is much heavier than the other. Such treatments do not capture 
genuine two-particle effects, like the two-particle relativistic corrections and hyperfine corrections. These corrections 
in general do not satisfy the scahng criterion and one expects the Crater property to be violated by these corrections. 
In passing, we note that the notion of Crater property can be generalized in a straightforward manner to any 
operator of the following form: 

0{R)^aS^{R-ri) + b5^{R^r2), (21) 

where a and 6 are arbitrary coefficients. This operator (i?) may correspond to physically interesting objects for 
particular choices of a and b; it is the charge density when (a, 6) = (61,62), the probability density of particle 
1 and 2 when (a, 6) — (1,0) and (0,1), respectively, and the mass density when (a, 6) = (7711,7712). Then the 
Hamiltonian or potential is said to exhibit the Crater property of charge/probability/mass distribution if and only if 
the charge/probability/mass distribution in the bound state of particle mass mi and 7712 is equivalent to the sum of 
the charge/probability/mass distributions produced in the limits ttii — > 00 and 7772 -^ cxd. As before, all these Crater 
properties are guaranteed by the same scaling criterion (S). 

In conclusion, we have provided a simple criterion to determine if the eigenfunctions of a given Hamiltonian have 
the Crater property. In particular, we have shown that neither the inclusions of fine structure corrections nor the 
switching from Schrodinger to Dirac formalism will destroy the Crater property. The author believes Crater must have 
foreseen the essential points of this paper — as in the conclusion of Ref. M he stated that "in general, the appearance 
of parameters in the potentials that are not dimensionless (in natural units) and do not depend on the reduced mass 
would not be of the correct type." As a consequence, this paper may be regarded as a concrete realization of this 
observation. 
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